Structural, electronic and magnetic properties of heterofullerene C48B12 
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Bonding, electric (hyper)polarizability, vibrational and magnetic properties of heterofullerene 
C48B12 are studied by first-principles calculations. Infrared- and Raman-active vibrational frequen- 
cies of C48B12 are assigned. Eight 13 C and 2 n B nuclear magnetic resonance (NMR) spectral signals 
of C48B12 are characterized. The average second-order hyperpolarizability of C48B12 is about 180% 
larger than that of Ceo- Our results suggest that C48B12 is a candidate for photonic and optical 
limiting applications because of the enhanced third-order optical nonlinearities. 
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I. INTRODUCTION 

In 1985, Kroto et al. [1] proposed the existence of C 6 o 
clusters in their graphite laser vaporization experiment. 
This proposal was subsequently confirmed in 1990 when 
Kratschmer et al. [2] reported a method for the mass 
production of Cgo in a carbon arc along with infra-red 
(IR) spectroscopic evidence for the Ceo carbon-structure. 
These pioneering works have stimulated extensive re- 
search into fullerenes [3,4], a new form of pure carbon, 
where an even number of three-coordinated sp 2 carbon 
atoms arrange themselves into 12 pentagonal faces and 
any number (> 1) of hexagonal faces. These carbon- 
cage molecules can crystallize into a variety of three- 
dimensional structures [2] and be doped in several dif- 
ferent ways [4]: endohedral doping, where the dopant 
is inside the fullerene cage; substitutional doping, where 
the dopant is on the fullerene cage; and cxohcdral doping, 
where the dopant is outside or between fullerene cages. 
It has been shown that doped fullerenes have remark- 
able structural, electronic, optical and magnetic proper- 
ties [4-6]. 

In 1995, the heterofullerene C59N was formed effi- 
ciently in the gas phase during fast atom bombardment 
mass spectroscopy of a cluster-opened N-methoxyethoxy 
methyl ketolactam [7]. The isolation and characteriza- 
tion of biazafullcrcnyl has opened a viable route for the 
preparation of C59N and other heterofullerenes in so- 
lution, leading to a number of detailed theoretical and 
experimental studies of C59N and heterofullerenes [4-6]. 
In 1991, the Smalley group [8] successfully synthesized 
boron-substituted fullerenes C6o- n B„ (1 < n < 6). Very 
recently, Hultman et al. [9] have successfully synthesized 
aza- fullerenes C6o-nN ra , formed by substituting carbon 
atoms in Ceo with more than one nitrogen atom, and the 
existence of a stable C48N12 aza-fullerene [9-12] was re- 
vealed. Stimulated by the high stability of C48N12, we 
have recently predicted that C48B12 [13] is also a stable 
heterofullerene and can be a promising component for 
molecular rectifiers, nanotube-based transistors and p-n 
junctions. 

In this letter, we further study the bonding, Mul- 



liken charges, electric (hyper)polarizability, vibrational 
and magnetic properties of C48B12. We characterize 13 C 
and 11 B NMR spectral lines of C48B12 and show how 
the boron-substitutional doping modifies the infrared and 
Raman spectra of the pristine Ceo- We also find that 
C48B12 exhibits enhanced second-order hyperpolarizabil- 
ity (enhanced third-order optical nonlinearity) and can 
compete with Cgo and aza-fullcrcnc C48N12 as a can- 
didate for photonic and optical limiting applications (for 
examples, data processing, eye and sensor protection, all- 
optical switching, and optical limiting) [6]. 

II. BONDING AND MULLIKEN CHARGE 

The geometry of C48B12, shown in Fig.l, was fully op- 
timized by using the Gaussian 98 program [14,15]. We 
have used the B3LYP [16] hybrid density functional the- 
ory (DFT) method, which includes a mixture of Hartree- 
Fock (exact) exchange, Slater local exchange [17], Bccke 
88 non-local exchange [18], the VWN III local exchange- 
correlation functional [19] and the LYP correlation func- 
tional [20], and a 6-31G(d) basis set. We consider of the 
form of C48B12 with such a dopant assignment: each pen- 
tagon has one boron atom and two boron atoms preferen- 
tially sit in one hexagon, while other forms are possible. 
The symmetry of C48B12 is found to be a Cj point group. 
The distances (or radii) R^ from the ith atom to the den- 
sity center of the molecule are listed in Table I. We find 
that C48B12, similar to C48N12 [11], is an ellipsoid struc- 
ture with 10 unique radii, while Geo h a s the same radius 
for each carbon atom (calculated Ri = 0.35502 nm, in 
excellent agreement with experiment [21]). 

The calculated net Mullikcn charges Q; of carbon and 
boron atoms in C48B12 are listed in Table I. Like C48N12 
[10,11], the heterofullerene C48B12 has two types of boron 
dopants in the structure: one with net Mulliken charges 
Qi = 0.1637 e and the other with Q ; = 0.1871 e. All 
carbon atoms in C48B12 have negative Q;. Although 
the Mulliken analysis can not predict exactly the atomic 
charges quantitatively, the sign of atomic charge can be 
estimated correctly [22]. From the Mulliken analysis, we 
see that boron atoms in C48B12 exist as electron donor 
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while carbon atoms act as electron acceptors. The calcu- 
lated charge of boron atom is consistent with the exper- 
imental result of the Smalley group [8] that an electron- 
deficient site was produced at the boron position on the 
cage. 




FIG.l: C48B12 geometric structure optimized with 
B3LYP/6-31G(d). The site numbers {5, 9, 14, 21, 26, 
30, 35, 39, 45, 50, 55, 60} are for boron atoms and the 
others for carbon atoms. 

Table I: B3LYP/6-31G(d) calculations of radius (R*, in 
nm) and net Mulliken charge (Qi, in e, where 1 e = 
1.6 x 10~ 19 C) for C 48 Bi2. 



Site Number {m} Atom Ri 



Qi 



{1, 13, 
{2, 12, 
{3, 11, 
{4, 15, 
{5, 14, 
{6, 18, 
{7, 19, 
{8, 20, 
{9, 21, 
{10, 17. 



16, 31, 

29, 32, 
28, 33, 
27, 34, 

30, 35, 

24, 42, 
23, 43, 
22, 44, 
26, 45, 

25, 41 



38, 51} 
37, 52} 
36, 53} 
40, 54} 

39, 55} 
48, 58} 
47, 57} 
46, 56} 
50, 60} 

, 49, 59} 



c 


0.35198 


-0.0385 


c 


0.34652 


-0.0079 


c 


0.35701 


-0.0036 


c 


0.37033 


-0.0712 


B 


0.36454 


0.1637 


C 


0.37116 


-0.0333 


C 


0.38035 


-0.0376 


C 


0.38252 


-0.0826 


B 


0.37953 


0.1871 


C 


0.37079 


-0.0761 



The optimized carbon-carbon (CC) and boron-carbon 
(BC) bond lengths in C48B12 are listed in Table II. We 
find that C48B12 has 6 unique BC bond lengths in the 
range of 0.15376 nm to 0.15899 nm and 9 unique CC bond 
lengths in the range of 0.13861 nm to 0.15014 nm. For 
Ceo molecule, the double C=C bond and single C-C bond 
lengths are 0.13949 nm and 0.14539 nm, respectively, in 
excellent agreement with experiment [21]. 

III. Electric (HYPER)POLARIZ ABILITY 

The static dipole polarizability (SDP) for hetero- 
fullerene C 48 Bi 2 is presented in Table III. The B3LYP 



results were obtained by using the Gaussian 98 program 
[14,15], while the LDA (local density approximation) re- 
sults were calculated by using the Amsterdam Density 
Functional (ADF) program [15,23,24]. The SDPs for 
C48N12 and Ceo listed in Table III are taken from our 
recent work [25]. For the B3LYP calculations, we use 
the valence-split basis set 6-31G(d) including the po- 
larization functions for boron and carbon atoms. The 
ADF program uses basis sets of Slater functions. In 
this work, we use a triple zeta valence plus polarization 
(TZP) augmented with field- induced polarization (FIP) 
functions of Zeiss et al. [26]. This basis set, TZP++ 
([6s4p2dlf] for carbon and boron atoms), was previously 
used for calculating the second-order hyperpolarizability 
7 of Ceo and its derivatives (for example, C58N2, C58B2 
and CssBN) and has been shown to produce reasonable 
(hyper)polarizabilities even with its small size [27]. Here, 
we only make a comparison between C48B12, C48N12 and 
Cgo- A comparison of Cgo with other theoretical and ex- 
perimental results can be found in review chapters [5,6] 
and recent work [25,27]. 

From Table III, we see that the LDA results are about 
20% larger than the corresponding B3LYP ones. This 
is expected since the basis set in the LDA calculation 
is larger and the LDA method, in general, predicts a 
larger polarizability than the B3LYP method [28]. Nev- 
ertheless, both B3LYP and LDA predict the same trends. 
The mean polarizability of C48B12 in the LDA (B3LYP) 
is about 12% ( 15%) larger than that of Cgo. 

The first hyperpolarizability /? of C48B 12 is zero due 
to inversion symmetry. The static second-order hyper- 
polarizability, 7, for C48B12, C48N12 and Cgo are pre- 
sented in Table IV. The 7 values for C4 8 N 12 and C 60 
listed in Table IV are taken from Ref . [25] , and the com- 
parison will only be made between these three molecules. 
For the calculations of the second-order hyperpolariz- 
ability, we use time-dependent DFT (TD-DFT) method 
as described in Ref. [27,29], i.e. finite- field differentia- 
tion of the analytically calculated first-order hyperpo- 
larizability. For all TD-DFT calculations, we used the 
RESPONSE code [15,30] implemented in the ADF pro- 
gram. [15,23,24]. 

We find for all components of the second-order hyper- 
polarizability for both C48B 12 and C4sN 12 a larger value 
than for Cgo- All 7 components for C48B12 are also larger 
than the corresponding components for C 4 sNi2 except for 
the j zzzz . This gives an average second-order hyperpo- 
larizability, 7 of C48B12, which is about 180% larger than 
that of Cgo- In contrast, the 7 value for C 4 sNi2 is about 
55 % larger than that of Cgo- The increase in the second- 
order hyperpolarizability is much larger for C48B12 than 
for C48N12 especially considering that there is no increase 
in volume. It has been experimentally shown that Cgo is 
a good optical limiter [6] because of its larger 7 value. 
Our present results imply that heterofullerene C48B12 
can compete with Cgo as an even better optical limiter 
because of its enhanced third-order optical nonlinearity. 
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Table II: B3LYP/6-31G(d) calculation of CC and BC bond lengths (in nm) in molecule C 48 Bi2. 



Bond Site Number Pairs for Bonding Bond Length 

~CC (1, 2) (12, 13) (16, 29) (31, 32) (37, 38) (51, 52) 0.14261 

BC (1, 5) (13, 14) (16, 30) (31, 35) (38, 39) (51, 55) 0.15376 

CC (1, 52) (2 , 31) (12, 38) (13, 29) (16, 37) (32, 51) 0.14078 

CC (2, 3) (11, 12) (28, 29) (32, 33) (36, 37) (52, 53) 0.15014 

CC (3, 4) (11, 15) (27, 28) (33, 34) (36, 40) (53, 54) 0.14871 

CC (3, 10) (11, 59) (17, 33) (25, 36) (28, 49) (41, 53) 0.13861 

BC (4, 5) (14, 15) (27, 30) (39, 40) (34, 35) (54, 55) 0.15667 

CC (4, 46) ( 8, 15) (20, 40) (22, 54) (27, 44) (34, 56) 0.13867 

BC (5, 42) ( 6, 35) (14, 48) (18, 55) (24, 30) (39, 58) 0.15576 

CC (6, 10) (17, 18) (24, 25) (58, 59) (41, 42) (48, 49) 0.14682 

CC (6, 7) (18, 19) (23, 24) (42, 43) (47, 48) (57, 58) 0.14080 

BC (7, 60) (9, 47) (19, 26) (21, 57) (23, 45) (43, 50) 0.15594 

CC (7, 8) (19, 20) (22, 23) (43, 44) (46, 47) (56, 57) 0.14570 

BC (8, 9) (20, 21) (22, 26) (44, 45) (46, 50) (56, 60) 0.15899 

BC (9, 10) (17, 21) (25, 26) (41, 45) (49, 50) (59, 60) 0.15563 

Table III: Static dipole polarizability (a, in nm 3 ) for C48B12, C48N12 and Cgo calculated with B3LYP/6-31G(d) and 
LDA/TZP++. The symmetry relation for Ceo gives a xx = ot vv — a zz and for C48B12 and C48N12 is a xx = a yy . 



B3LYP/6-31G(d) LDA/TZP++ 
Molecule a xx a zz a xx a zz Ref. 

C^ 0.0695 0.0695 0.0847 0.0847 [25J 

C 48 Ni2 0.0666 0.0675 0.0793 0.0815 [25] 

C48B12 0.0822 0.0804 0.0958 0.0939 this work 

Table IV: The static second-order hypcrpolarizabilities (7, in a.u., with 1 a.u. = 6.235378 xl0~ 65 C 4 m 4 J- 3 ) for 
C48B 12 , C4 8 N 12 and C 60 calculated by using LDA and a TZP++ basis set. The average second-order hyperpolar- 
izability is given by 7 = ^ J2i.j(lan + lijij + lijji)- Tne symmetry relations of the molecule gives ~j xxxx = 'Jyyyy, 
"f xxzz — r )yy zz and ^y zzxx — r y ZZ yy- 



Molecule "/ xxxx r ) xxyy ^i zzzz r ) xxzz r ) zzxx 7 Ref. 
C 60 137950 45983 137950 45983 45983 137950 [25] 
C 48 Ni2 188780 62880 232970 85120 84790 215222 [25] 
C 48 Bi2 470190 156840 214300 116800 118090 387628 this work 



IV. IR AND RAMAN SPECTRA 

Using the Gausian 98 program [14,15], we first opti- 
mize the geometry of C48B12 and Ceo with the B3LYP 
method and 3-2 1G basis set. Then, we calculate the vi- 
brational frequencies of C48B12 and Ceo with the same 
method and basis set. Our results for C@o arc in agree- 
ment with experiment [31,32]. Ceo has totally 46 vibra- 
tional modes [4]. Since C48B12 has lower symmetry (d) 
than Ceo, we find 174 independent vibrational modes for 
C48B12: 87 non-degenerate IR-active modes with a u sym- 
metry and 87 non-degenerate Raman-active modes with 
a g symmetry. 

We also calculate IR intensities Im and Raman scat- 
tering activities O rama „ at the corresponding vibrational 
frequencies for both Cgo and heterofullerene C48B12. The 
results are shown in Fig. 2 and Fig. 3. Since experimental 
IR and Raman spectroscopic data do not directly indi- 
cate the specific type of nuclear motion producing each 



spectroscopic peak, we do not give here the normal mode 
displacement for the vibrational frequencies. 

For Ceo, we note that its IR spectrum is very sim- 
ple. Namely, it is composed of 4 IR spectroscopic sig- 
nals with t\ u symmetry. The IR intensities for Cgo cal- 
culated with B3LYP/3-2G agree reasonably with the in 
situ high-resolution FTIR spectrum of a Ceo fihn mea- 
sured by Onoe and Takeuchi [33]. However, the IR spec- 
trum of C48B12 is not so simple, exhibiting 87 IR spectro- 
scopic signals, with the stronger IR spectroscopic signals 
mainly in the high-frequency region. Here we discuss 
several vibrational bands in the IR intensity of C48B12 
to be compared with future experimental identification: 
(i) the strongest IR spectroscopic signal with IR inten- 
sity of 229912 m/mol is determined at the high frequency 
1356 cm -1 ; (ii) three strong modes are observed at fre- 
quencies v — 783 cm -1 , 1031 cm -1 and 1546 cm" 1 with 
IR intensities of 85828 m/mol, 119322 m/mol and 114878 
m/mol, respectively; (iii) five intermediate modes appear 
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at frequencies v = 730 cm -1 , 1018 cm -1 , 1106 cm -1 , 
1217 cm -1 and 1303 cm -1 with IR intensities of 58303 
m/mol, 69950 m/mol, 55094 m/mol, 69160 m/mol and 
53856 m/mol, respectively. 
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Fig.2: B3LYP/3-21G calculation of IR-active vibra- 
tional frequencies (u, in cm -1 ) and IR intensities (Iir, 
in 10 3 m/mole) of C48B12 (open circles) and Cqo (solid 
lines). 
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Fig. 3: B3LYP/3-21G calculations of Raman-active 
frequencies {y, in cm -1 ) and Raman scattering activi- 



ties ( n r 



in 10- 14 m 4 /kg) of C 48 Bi 2 . The solid and 



dot-dashed lines are the unpolarized and polarized Ra- 
man spectral lines of Cgo , respectively. Filled circles and 
open squares are non-degenerate unpolarized and polar- 
ized Raman-active modes, respectively. 

As shown in Fig. 3, Ceo has two non-degenerate polar- 
ized Raman spectroscopic signals with a g symmetry and 
8 fivefold-degenerate unpolarized ones with h g symme- 
try. The strongest Raman spectroscopic signals in Ceo 
are the two a g modes. The calculated results for Cqo are 
in excellent agreement with experiment [32]. In contrast, 
for C48B12, we observe 10 non-degenerate unpolarized 
and 77 non-degenerate polarized Raman spectroscopic 
signals with a g symmetry. The Raman spectrum sepa- 
rates into high-frequency (1000 cm -1 to 1700 cm -1 ) and 
low- frequency (200 cm -1 to 900 cm -1 ) regions, similar to 
those of Ceo- The strong Raman spectroscopic signals in 
C48B12 are the non-degenerate polarized modes. Three 



strong Raman bands of C48B12 are: (i) The strongest 
Raman spectroscopic signal in C48B12 appears at high 
frequency v = 1394 cm -1 with Q r aman — 5035 x 10~ 14 
m 4 /kg; (ii) four strong Raman spectroscopic signals are 
located at v = 1014 cm -1 , 1053 cm -1 , 1353 cm -1 
and 1482 cm" 1 with fl raman = 1130 x 10~ 14 m 4 /kg, 
2094x 10~ 14 m 4 /kg, 1568x 10~ 14 m 4 /kg and 2149x 10~ 14 
m 4 /kg, respectively; (hi) five intermediate strong Ra- 
man spectroscopic signals are observed at v = 1105 
cm -1 , 1166 cm -1 , 1214 cm -1 , 1336 cm" 1 and 1518 cm -1 
with n raman = 618 x 10~ 14 m 4 /kg, 644 x 10~ 14 m 4 /kg, 
590xl0~ 14 m 4 /kg, 601 xl0~ 14 m 4 /kg, 537xl0~ 14 m 4 /kg 
and 887 x 10~ 14 m 4 /kg, respectively. 

V. SECOND-ORDER MAGNETIC RESPONSE 

There are a number of theoretical methods for cal- 
culating the second-order magnetic response proper- 
ties of molecules. In this paper, we use both the 
gauge-including-atomic-orbital (GIAO) method and the 
continuous-set-of-gauge-transformation (CSGT) proce- 
dure [34], which is implemented in the Gaussian 98 pro- 
gram [14,15], to predict the NMR shielding tensors a 
of C48B12. In high- resolution NMR, the isotropic part 
(Ti SO of a is measured by taking the average of a with 
respect to the orientation to the magnetic field, i.e., 
{a xx + v yy + cr zz )/3, where and a zz 

are the principal axis values of a. The results calcu- 
lated by using B3LYP hybrid DFT and restricted Hartree 
Fock (RHF) theory are summarized in Table V. We find 
that C 48 Bi2 has 8 13 C and 2 n B NMR spectral signals, 
indicative of the 10 unique sites in the C48B12 struc- 
ture. In contrast, Ceo has only one 13 C NMR spec- 
tral signal, for example, with Ui SO — 50.5 ppm (parts 
per million) and 54.7 ppm obtained with B3LYP/6- 
31G(d):GIAO and RHF/6-31G(d):GIAO, respectively. 



The 13 C NMR chemical shift {5 = a. 



TMS 



sample 



) with 



respect to the reference tetramethysilane (TMS) for Ceo 
is, for example, S = 133.3(135.7) ppm for B3LYP/6- 
31G(d):GIAO (B3LYP/6-31G(d):CSGT), about 9 (7) 
ppm difference from experiment (5 = 142.7 ppm [35]), 
but S = 140.4(141.7) ppm for RHF/6-31G(d):GIAO 
(RHF/6-31G(d):CSGT) which is in good agreement with 
experiment. The results for Cqo show that the DFT 
method does not provide systematically better NMR re- 
sults than RHF. This is due to the fact that no current 
functionals include a magnetic field dependence [34] . For 
Ceo, the CSGT procedure provides better NMR results 
than the GIAO procedure, but takes more CPU time (see 
Table V) than the GIAO procedure when compared to 
experiment. 

VI. SUMMARY 

In summary, we have performed first-principles calcu- 
lations of bonding, Mullikcn charges, dipole polarizabil- 
ity, hypcrpolarizability, vibrational frequencies, IR in- 
tensities, Raman scattering activities and second-order 
magnetic response properties of hctcrofullcrcnc C48B12. 
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Eighty-seven independent IR-active and 87 independent hibits enhanced third-order optical nonlinearity, which 
Raman-active vibrational modes for C48B12 are assigned. implies potential applications of C48B12 in photonics and 
Eight 13 C and two n B NMR spectral lines for C48B12 are optical limiting, 
characterized. Compared to C60 and C48N12, C48B12 ex- 
Table V: B3LYP/6-31G(d) and RHF/6-31G(d) calculations of the absolute isotropy, Ui SO in ppm (parts per million), 
of the nuclear magnetic shielding tensor a for atoms in C48B12, Ceo and tetramethysilane (TMS) found by using both 
GIAO and CSGT methods. The CPU times (in hours) for Ceo cases are shown. 
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